Scanning tunneling spectroscopy of a magnetic atom on graphene in the Kondo regime 
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The Kondo effect in the system consisting of a magnetic adatom on the graphene is studied. By 
using the non-equilibrium Green function method with the slave-boson mean field approximation, 
the local density of state (LDOS) and the conductance are calculated. For a doped graphene, the 
Kondo phase is present at all time. Surprisingly, two kinds of Kondo regimes are revealed. But for 
the undoped graphene, the Kondo phase only exists if the adatom's energy level is beyond a critical 
value. The conductance is similar to the LDOS, thus, the Kondo peak in the LDOS can be observed 
with the scanning tunneling spectroscopy. In addition, in the presence of a direct coupling between 
the STM tip and the graphene, the conductance may be dramatically enhanced, depending on the 
coupling site. 
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Graphene, a single-layer hexagonal lattice of carbon 
atoms, has recently attracted a great deal of attention 
due to its unique properties and potential applications.^ 
The graphene has an unique band structure with a 
linear dispersion near the Dirac-points, such that its 
quasi-particles obey the two-dimensional Dirac equation 
and have the relativistic-like behaviors. For a neutral 
graphene, the Fermi level passes through the Dirac points 
and the density of states (DOS) at the Fermi face van- 
ishes, so that its transport properties are greatly deviated 
from that of a normal metal. 

The Kondo effect has been paid great attention in con- 
densed matter community over the yearsj^^ It is a pro- 
totypical many-body correlation effect involving the in- 
teraction between a localized spin and free electrons. The 
Kondo effect occurs in the system of a magnetic impu- 
rity embedded in a meta&^ or a quantum dot coupled 
to the leadsj£i£i£ in which the magnetic impurity or the 
quantum dot acts as a localized spin. At low tempera- 
ture, the localized spin is screened by the free electrons, 
resulting a spin singlet and a very narrow Kondo peak 
located near the Fermi surface in the local DOS (LDOS). 

For a conventional metal, it has a finite DOS at the 
Fermi surface, thus, at zero temperature the Kondo ef- 
fect can occur with any weak exchange interaction J ex 
between the magnetic impurity and the free electrons. 
But in certain special systems, such as the nodal d-wave 
superconductors or the Luttinger liquids, the DOS van- 
ishes at the Fermin surface with power law behavior. In 
these systems, there exists a critical exchange coupling 
J c , and the Konde effect only occurs when the exchange 
interaction J ex > J c £ 

The graphene, which was successfully fabricated in re- 
cent years, has an unique band structure pi For the un- 
doped neutral graphene, its DOS is directly proportional 
to the energy |e| and the DOS vanishes at the Fermi sur- 
face, so it has an unconventional Kondo effect. Very 
recently, some studies have investigated the Kondo ef- 
fect in the graphene and a finite critical Kondo coupling 



strength was revealedi 10 ! 11 By varying the gate voltage, 
the charge carriers of graphene can be easily tuned exper- 
imentally. Then the Fermi level can be above or below 
the Dirac points, and the DOS depends on the energy e 
in the form |e — eo|, with the Dirac-point energy denoted 
by £q. In this situation, the DOS is quite small but fi- 
nite at the Fermi surface. In addition, there exists a zero 
point in DOS that is very close to the Fermi surface. How 
is the Kondo effect affected by this unique energy band 
structure? 




FIG. 1: (Color online) The schematic diagram for a magnetic 
adatom on the graphene probed by a STM tip. 

In this Letter, we study the Kondo effect in the 
graphene. We consider the model of a magnetic adatom 
on the graphene and to investigate its scanning tunneling 
spectroscopy by attaching a STM tip near the magnetic 
adatom (see Fig.l). Taking the slave-boson (SB) mean 
field approximation and using the non-equilibrium Green 
function method, the expressions of the LDOS and con- 
ductance are obtained. Our results exhibit that for the 
doped graphene the Kondo phase always exists, however, 
two kinds of Kondo regimes emerge. For the undoped 
neutral graphene, the Kondo phase only exists when the 
energy level td of the adatom is larger than a critical 
value e dc . 

The above mentioned system can be described by the 
Anderson Hamiltonian. Here the magnetic atom is mod- 
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eled by a single level td with the spin index a =1,1 and 
a constant on-site Coulomb interaction U. In the limit 
U — > oo, the double occupancy is forbidden. Then the 
SB representation can be applied, and the Hamiltonian of 
the Anderson model describing the system is transformed 
to the following form^ 

H = €dha + toa\ a a la + t Y a \* a ]<? 

<J i,a {ij),tr 

+ Y £fcC L c ^ + Y yLc\jab ] + t LR c\ a a Ja + ff.C.J 

k,a k,<T 

+ fafP>*o° + H.c.) + \(fitb + - 1) (1) 



where n CT = ai CT and Cka are the annihilation op- 

erators at the discrete site i of the graphene and the tip 
of STM, respectively. Here the graphene is assumed to 
be two-dimensional and without disorder. In Hamilto- 
nian (1), the graphene is described by the tight-binding 
model and only the nearest neighbor hopping t is consid- 
ered. f a and b are the pseudofermion operator and SB 
operator of the magnetic atom. The last term represents 
the single-occupancy constraint b^b+^2 a h a — 1 with La- 
grange multiplier A. The magnetic atom couples to the 
graphene at the lattice site with the coupling coefficient 
tR. The STM tip couples to both of the magnetic atom 
and graphene, and and tLR are their coupling coeffi- 
cients. Here we consider that the STM tip only contacts 
to a lattice site J of the graphene. Following, we take 
the standard SB mean field approximation in which the 
operator b is replaced by the c-number (b) . This approx- 
imation is a qualitative correct for describing the Kondo 
regime at zero temperature^ For convenience, we also 
introduce the renormalized energy id = td + A and hop- 
ping elements =ti, (b) and tR = tR (b). 

Next, we consider that a bias V is applied between 
the STM tip and the graphene to obtain the current I 
flowing from STM to the graphene and the LDOS around 
the magnetic atom. By using the non-equilibrium Green 
function method, the current I and LDOS are^ 
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t L Gf L (u)+t LR G< L (uj) 
LDOSH = -(l/7r)ImG^(w), 



(2) 
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where G < and G r are the standard lesser and retarded 
Green's functions, and G^^ L , JL in Eq.(2) and (3) are 

the elements of the matrix G < / r . Under the SB mean 
field approximation, G*" and G r can be solved from the 
Keldysh and Dyson equations: G< = (1 + G r E r )g<(l + 
S°G a ) and G r = (g r - £ r )-\ where £ r ' a are the 
self-energies. T, r Lf = Y7 fL = t L , ££ = W Qf = t R , 
j = S T 7L = tLR, and other elements of S r are zero. 
g r / < are the Green's functions of the isolated STM tip, 
magnetic atom, and graphene (i.e. while tr, = t R — 
= 0). They can easily be obtained: g r LL (uj) = 
Efeffc = -inPL, 9 r f f(u) = l/(w -e d + iO + ), g^ L {u) = 



sLfe„( w ) 



(uj - e ) 2 - \(f)\ 2 \ 4>* ^ - eo 



(5) 



where p L is the DOS of the STM tip. f L (w) and 
fniu}) are the Fermi-Dirac distribution functions in the 
isolated STM tip and graphene, respectively. gy < (w) 
are the Green's functions of the isolated graphene, and 
i = [ix,iy) and j = (jx,jy) are the indices of the unit cell. 
Since there are two carbon atoms per unit cell, gy is a 
2x2 matrix. In Eq.(4) and (5), a is the lattice constant 
of the graphene, the integral dk x dk y runs over the first 



Brillouin zone, and <f> = t(l + e * + e 



-iak x +iak y 



)■ 



At end we need to self-consistently calculate two un- 
knowns (b) and A with the self-consistent equations: 



(fe) 2 +n T +n i = l 



(6) 



27m 



du> 



^ (t L G< f (oj) + t R G< f (w)) + X(b) 2 = 0(7) 



where n a — f(doj/2ir)ImGjf(uj) is the electron occupa- 
tion number in the spin state a of the magnetic atom. 
In the numerical calculations, we take the hopping en- 
ergy t = 1 as the energy unit. The chemical potential 
Hr of the graphene is set to be zero as the energy zero 
point, then the bias voltage V = \xl- The temperature 
T is assumed to be zero since ksT is much smaller than 
t = 2.75eV in the real experiments. 

First, we assume the system to be decoupled to the 
STM tip (i.e. at tL = tLR = 0) and focus on the Kondo 
effect. Fig. 2 shows magnetic-atom LDOS for different e<2 
and eo- Here the Dirac-point energy e can be experimen- 
tally tuned by the gate voltage. For the magnetic atom, 
the level is fixed. But for an artificial atom, such as 
a quantum dot, is also tunable.- As shown in Fig. 2a, 
there exists three peaks in the LDOS. Two peaks locate 
at u> = ±t, due to the peaks in the DOS of the graphene. 
The third peak is the Kondo peak, which is always lo- 
cated around hr — regardless of the parameters of e^, 
eo, and tR. Fig. 2b magnifies the curves of Fig. 2a in the 
vicinity of u> = 0. It clearly shows that the Kondo peak 
is much higher than the other two peaks at u> = it. The 
lower the level td is, the higher and sharper is the Kondo 
peak and the closer it is to (1r = 0. In particular, for 
the undoped neutral graphene with e = 0, the Kondo 
peak is greatly deflected from (see Fig. 2b) although it 
seems to be near on a large scale (see Fig. 2a). The 
LDOS is zero at u> = since the graphene's DOS van- 
ishes at this point. This characteristics is much different 
from the conventional Kondo effect, in which the LDOS 
is quite large at u) = 0. 

Next, the energy eo is tuned to be away from zero. 
While eo < 0, the carriers in graphene are electron-like, 
but they are hole- like for e > 0. With e decreasing 
from 0, the Kondo peak is gradually getting closer to 
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FIG. 2: (Color online) LDOS vs. the energy uj with tz, = 
thR = and tR = 0.9t. In (a), the Dirac-point energy eo = 0. 
(b) is the enlarged view of the Kondo peaks in (a), (c) and 
(d) are at different eo for the level e^ = — 0.6i. 



FIG. 3: ( Color online) (a) Tk vs. the level td and (b) e^ vs. 
ed at tR — 0.9t for different energy eo. (c) Tk as functions of 
ed for eo = at different coupling tR. (d) Phase diagram of 
the eo = in the parameter space of (ed, tR). 



and is slightly widened (see Fig. 2c). It recovers the 
shape of conventional Kondo peak at large negative eo 
(e.g. e = -O.Oli or -0.02t in Fig.2d). On the other 
hand, with eo increasing from 0, the Kondo peak departs 
gradually from and its shape never recovers to that of 
the conventional Kondo peak even at large eo. 

From the position w\ and half-width Tr of the Kondo 
peak, we can estimate the Kondo temperature byi T k = 

*J 'ujf + f 2 R (cji), where the peak position uj\ is obtained 
from the equation uj — id — Re [t R gQ (uj)] = and the 
peak half-width Tr — — 2^Im[<7o (w)]. Fig. 3a shows the 
Kondo temperature Tk versus the level ed at a fixed tR 
for different energy eo- Tk possesses the following char- 
acteristics: (i) Tk decreases monotonically with reducing 
of ed- When ed is approaching — oo, Tk goes to zero, (ii) 
For eo 7^ 0, there exists two different Kondo regimes. 
When ed is larger than a critical value edc, Tk is almost 
linearly dependent on e^ — edc, so that Tk drops rapidly 
while ed is near edc (see Fig. 3a). This characteristics is 
very different comparing to the conventional Kondo ef- 
fect, and it origins from the linear DOS of the graphcnc 
and vanish-ness of the DOS at the Dirac point. On the 
other hand, when ed < edc, Tk oc e £d (i.e. lnTft oc ed as 
shown in Fig. 3a). In this case the system's behavior is 
similar to that of the conventional Kondo effect. Since 
ksTk <C |eo — ijlr\ in this regime, the DOS at the Fermi 
surface is approximatively constant on the energy scale 
ksTk- So it recovers the conventional Kondo behavior, 
(iii) The behavior of Tk for eo < (i.e. the electron doped 
graphene) and eo > (i.e. the hole doped graphene) is 
similar. 

Fig. 3b shows the renormalizcd level id which deter- 
mines the position of the Kondo peak. While eo < 0, the 
renormalized level id is always positive and it decreases 
monotonically with reducing of ed- id approaches zero as 



ed — ► — oo. On the other hand, while eo > 0, id can be 
negative. The curve of id~Sd has a negative minimum. 
But id — * is still true as e^ — > — oo. 

Now we focus on the case of eo =0. When the Dirac- 
point energy eo approaches to 0, the zero-DOS point ap- 
proaches to the Fermi surface, then both of the Kondo 
temperature Tk and the renormalized level id goto zero 
in the whole region of e^ < edc- At eo = (i.e. for 
the undoped neutral graphene), Tk and id reach zero at 
ed = edc- On the ed < edc side, the self-consistent equa- 
tions (6) and (7) do not have a solution. In other words, 
the system cannot be in the Kondo phase for ed < edc- 
On the other side with e^ > e^, the Kondo phase does 
exist. The Kondo temperature Tk is roughly proportional 
to ed — edc (see Fig. 3c). The critical value edc depends 
on the coupling strength tR between the magnetic atom 
and the graphene. The larger tR is, the larger |e<j c | is (see 
Fig. 3c). From Tk = or e~d — at e^ = ed c and with the 
aid of Eq.(7), the critical value edc can be analytically 
obtained: 

7T J-oo U 

The above integral over uj is about 1.408, hence e<j c w 
0.896i|j,. Fig. 3d shows the phase diagram in the param- 
eter space of (ed, £r)- The curve gives the critical value 
edc versus the coupling strength tR. In the top right re- 
gion, the Kondo phase emerges at low temperature. But 
in the bottom left region, the system cannot enter into 
the Kondo state even at zero temperature. 

Finally, we investigate the conductance G (G = 
dl/dV) while the STM tip is coupled. The coupling be- 
tween the STM tip and the adatom is normally much 
weak compared to that between the graphene to the 
adatom, so we set t^ and t^R on the order of O.Olt. Fig. 4 
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FIG. 4: (Color online) Conductance G vs. the bias V for 
different energy eo with tn — 0.9t, tz, = 0.002t, and td = 
—O.Gt. In (a), tLR = 0, and in (b), (c), and (d), tLR = 0.03t. 
The contact point J of the STM tip and the graphene in (b), 
(c), and (b) are 0, 1, and 2, respectively. 

shows the conductance G versus the bias V in the Kondo 
regime. While without the direct coupling (tLR = 0), 
the curve of G-V is very similar to the curve of LDOS-w 
(see Fig. 4a and 2c) since for the weak-coupled STM tip, 
the conductance G is proportional to the LDOS of the 
adatom. 14 When the direct coupling is present (tLR ^ 0), 
the coherence between the direct path and the path pass- 
ing the magnetic adatom occurs, which usually leads to 
the Fano resonance J£ But in the present graphene sys- 
tem, the Fano resonance does not appear anywhere, the 
reason is that the transmission coefficient of the direct 
path also depends on the energy of the incident electron. 
In Fig.4b, 4c, and 4d, J of the contact point of the STM 
tip to the graphene is (same contact point with the 
adatom), 1 (the nearest neighbor site), and 2 (the next 



nearest neighbor site), respectively. For the same and the 
next nearest neighbor contact points, the conductance G 
is almost unaffected by the opening of the direction path, 
although tLR = 0.03t is much larger than tL = 0.002£. 
Because that in the Kondo regime, a high Kondo peak 
emerges in the LDOS of the adatom, so the path passing 
the adatom is large and plays the dominant role. How- 
ever, when the contact point is at the nearest neighbor 
site, the conductance is greatly enhanced by the open- 
ing of the direct path (see Fig. 4c, note the different scale 
comparing with the other plots in Fig. 4) since the Kondo 
singlet state is quite extensive. For the graphene system, 
the sharp Kondo peak not only appears in the LDOS of 
the adatom, but also in the LDOS of all odd-neighboring 
sites. The large LDOS in the odd-neighboring sites and 
tLR > tL lead that the direct path is dominant, and the 
conductance G is enhanced. 

In summary, the Kondo effect in the system of a mag- 
netic adatom on the graphene and its scanning tunnel- 
ing spectroscopy are theoretically studied. For a doped 
graphene, the Kondo phase always exists at zero tem- 
perature regardless of n-type or p-type doping, and it 
exhibits two Kondo regimes: the conventional and un- 
conventional. For the undoped neutral graphene, the 
Kondo phase only exists while the adatom's energy level 
is above a critical value. These Kondo characteristics 
can be observed from scanning tunneling spectroscopy 
since the conductance versus the bias is similar to the 
LDOS of the adatom. In addition, if a direct coupling 
between the STM tip and the graphene is present, the 
conductance may be greatly enhanced, depending on the 
coupling site. 
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